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ABSTRACT 

We calculate the structure of a relativistic shock wave in which the internal 
energy of the shocked fluid is radiated away on a time scale much shorter 
than the characteristic shock propagation time. The shock is assumed to 
move through a uniformly magnetized, neutral plasma consisting of protons, 
electrons, and positrons, and allowance is made for the possible production of 
electron/positron pairs in the shock itself. The radiation mechanism is taken 
to be synchrotron and inverse- Compton emission (the latter involving both 
synchrotron-produced and externally supplied seed photons) by the electrons 
and positrons. We simplify the discussion by considering a shock in which the 
magnetic field is transverse to the direction of propagation and focus attention 
on the properties of the radiative zone that forms behind the shock transition. 
In particular, we investigate the possibility that the compression induced by 
the cooling of the gas amplifies the magnetic pressure until it reaches (and 
ultimately exceeds) equipartition with the thermal pressure (which, in turn, 
limits the overall compression). We show that, if a significant fraction of the 
postshock thermal energy is deposited in the electron/positron component, then 
a considerable portion of the emitted radiation will come from regions of strong 
magnetic field even if the field immediately behind the shock transition is well 
below equipartition. This field amplification mechanism may be relevant to the 
production of synchrotron flares in blazars, miniquasars, and gamma-ray burst 
sources. We consider the latter application in some detail and show that this 
process may play a role in the prompt 7-ray and possibly also the optical "flash" 
and radio "flare" emission, but probably not in the afterglow. 

Subject headings: galaxies: active — gamma rays: bursts — radiation mechanisms: 
nonthermal — shock waves 
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1. Introduction 

Propagating relativistic sliock waves are the leading candidates for the origin of the 
emission from gamma-ray bursts (GRBs) and their afterglows (e.g., Piran 1999, 2000; 
Meszaros 2000) as well as of the spectral flares and apparent superluminal motions exhibited 
by the blazar class of active galactic nuclei (AGNs; e.g., Ulrich, Maraschi, & Urry 1997) and 
their "miniquasar" Galactic counterparts (e.g., Mirabel & Rodriguez 1999; [Kaiser, Sunyaev, 



Spruit 200C ). The shocks are thought to be associated with relativistic, jet-like outflows 



that are driven from the vicinity of a compact object: a massive black hole in blazars, a 
stellar-mass black hole in miniquasars, and a stellar-mass black hole or a rapidly rotating 
neutron star in GRBs. Q In all cases except GRB afterglows, the shocks are believed to be 
internal, possibly resulting from the collision of "shells" (the products of nonsteady ejection 
episodes) inside the jet (e.g., Rees 1978; [benimore et al. 1996| ; |Sari fc Piran 199TBD . In the 



case of GRB afterglows, the most common interpretation is in terms of shocks that form at 
the interface between the relativistically outflowing material and the surrounding medium 
(with the bulk of the observed emission arising in the "forward" shock that propagates 
into the ambient medium, whereas the "reverse" shock that is driven into the ejecta gives 
rise to an early "optical flash" and associated "radio flare"). The radiation in all three 
types of sources is inferred to be nonthermal, with the dominant emission mechanisms 
most commonly invoked being synchrotron and inverse Compton. In particular, both the 
7-ray burst emission and the radio-through-X-ray afterglow radiation in GRBs, as well 
as the radio-through-X-ray emission in blazars and the radio-through-optical radiation 
in microquasars, has been attributed to synchrotron radiation (or, in certain models of 
GRBs, a closely related mechanism; e.g., [Smolsky fc Usov 2000] , [Medvedev 2000] , P^yutikov 



|fc Blackman 2001| ). In the case of blazars, the measurement of high 10%) optical linear 



polarization in some objects directly supports the synchrotron interpretation (e.g.. Angel fc 
Stockman 1980). In the case of GRBs there are already two confirmed detections of optical 
polarization (Covino et al. 1999; Wijers et al. 1999; Rol et al. 2000), which, despite the 
relatively low measured values (~ 1 — 3%), appear to be consistent with a synchrotron 
origin. In both blazars and microquasars the case for synchrotron radiation is supported 
also by measurements of radio polarization. 

Synchrotron radiation requires the presence of charged particles (typically electrons or 
electron-positron pairs) with relativistic random velocities that move in a magnetic field. 
The magnitude of the energy densities of the random relativistic electron motions and of the 



■^In at least some GRBs there are, however, indications that the outflow may have a fairly large opening 
angle or even be spherical. 
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magnetic field are parametrized in terms of tlieir ratios to tlie total thermal energy density 
in the emission region by ee and e^, respectively. In most cases, the values of these two 
parameters cannot be accurately determined from the observational data. However, in the 
case of blazars there are indications that the emitting gas is in approximate equipartition 
between the particle and magnetic energy densities (i.e., eg ~ — e.g., Readhead 1994; 
Bower & Backer 1998). This conclusion is also consistent with models of the nonthermal 
flares in microquasars (e.g.. Kaiser et al. 2000). In the case of GRBs, if the observed 7-ray 
emission indeed represents synchrotron radiation, then considerations of radiative efficiency 
(based on estimates of source energetics) imply that eg, and likely also e^, are not much 
smaller than 1 (e.g., ^ari fc Piran 199"7a| ; Piran 1999). In the case of GRB afterglows, there 



have been a few instances in which data were available for deriving more direct estimates 
of these parameters. In particular, theoretical fitting of the broadband spectrum of the 
afterglow of GRB 970508 has yielded ~ 0.1 - 0.6 and €3 ~ 0.01 - 0.1 ( |Wijers fc Galaim 
TOTt Mranot, Piran, fc Sari 1999| ; [Chevalier fc Li 2000| ). 



The inference that ee and are not much smaller than 1 in GRBs and at least 
some afterglows has important implications to the nature of the emitting regions in these 
sources. A comparatively large value of implies a high efficiency of shock acceleration of 
relativistic electrons in an electron-proton plasma or else a high electron-positron (e^) pair 
content in the postshock region (reflecting the composition of the preshock gas or arising 
from copious pair production in the shock). A comparatively large value of implies either 
a highly magnetized preshock medium or else an efficient mechanism of field amplification 
in (or behind) the shock. If eg ~ 1 and the radiative cooling time is much shorter than 
the dynamical timeQ (which, for synchrotron radiation, would require also not to be too 
small), then the shock will be radiative: it will contain a radiative cooling layer, where 
most of the internal energy is radiated away, that will directly infiuence its dynamics. 
Nonrelativistic, thermally emitting shocks of this type are prevalent in the interstellar 
medium (e.g., praine fc McKee 1993| ), but their properties in relativistic, nonthermal 



sources have not yet received much attention in the literature. 

A radiative shock wave generally consists of a "shock transition," where the bulk of 
the kinetic energy of the preshock gas is dissipated, and a "radiative zone," where most 
of the dissipated energy is radiated away. The energy dissipation scale is typically much 
smaller than all other relevant length scales, so the shock transition can be treated as an 
infinitesimal front at the upstream end of the shock wave. The radiative zone, on the other 



*The dynamical time is the time for a significant expansion of tlie emitting region, and is thus also the 
characteristic time for adiabatic energy losses. It is typically of the order of the shock distance from the 
origin divided by the instantaneous shock speed. 
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hand, has a finite width, corresponding to the coohng length over which the postshock 
gas radiates away a significant fraction of the dissipated energy. After passing through 
the shock transition, the gas is usually compressed by a factor of a few (for Newtonian or 
mildly relativistic shocks), and in the absence of significant radiation this is the density 
enhancement that also characterizes the far-downstream (or postshock) region. However, 
it is a well-established result in the theory of interstellar shocks that the asymptotic 
compression can be far larger if the shock is radiative. This is a consequence of the fact that 
the total (thermal plus magnetic) pressure behind the shock transition remains roughly 
constant, so that, as long as magnetic effects are not important, the particle density n 
increases inversely with decreasing temperature T. 

A related important result for radiative interstellar shocks is that, if the gas has a 
"frozen in" magnetic field B with a significant component {B±) perpendicular to the shock 
normal, then magnetic stresses become progressively more important as T decreases behind 
the shock, a reflection of the fact that the magnetic pressure, associated with B±, scales 
as n^. If the preshock magnetic field is sufficiently strong, then the magnetic pressure in 
the cooling layer will reach equipartition with the thermal pressure at some point within 
the radiative zone. Beyond this point the gas will undergo little additional compression, 
and the cooling will proceed at nearly constant density rather than at constant thermal 
pressure (e.g., HoUenbach & McKee 1979). Since the cooling time is proportional to 1/n (a 
consequence of the fact that the cooling function A, the energy radiated per unit volume 
and per unit time, scales as for a thermally emitting gas), the inhibition of further 
compression by the magnetic stress results in a lower cooling efficiency in this part of the 
radiative zone. 

One can anticipate that the basic properties of radiative hydromagnetic shocks in 
the interstellar medium (ISM) would carry over to the relativistic, nonthermal regime. In 
particular, one can expect a significant enhancement of the compression as well as of the 
downstream magnetic field in a radiative shock. The two cases are, however, expected 
to differ in detail. In particular, the cooling function of a synchrotron-emitting shocked 
gas satisfies A oc B^rif. oc nl (assuming a predominantly transverse magnetic field behind 
the shock), where rie is the number density of electrons (and positrons). The explicit 
density dependence of A is thus stronger than in a thermally emitting gas, which suggests 
that the dynamical infiuence of the magnetic field on the structure of the shock would be 
even more pronounced in the nonthermal case. This could have important implications 
to astrophysical sources in which such shocks occur. For one thing, it may provide a 
mechanism for amplifying the magnetic field up to equipartition strengths in cases where 
it is initially comparatively weak. Furthermore, as the processes of field amplification and 
synchrotron emission are mutually reinforcing in the radiative zone, this mechanism may 
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naturally lead to synchrotron emission regions that are characterized by a high value of e^. 
There is even the possibility, if inverse- Compt on cooling also plays a role, that the bulk 
of the emitted synchrotron radiation would come from regions where the magnetic field is 
already at (or near) equipartition.^ 

In this paper we examine the aforementioned effects in a quantitative manner by 
studying the properties of relativistic, hydro magnetic, nonthermally emitting shocks. In §0 
we describe the framework of our analysis and calculate the shock jump conditions. The 
structure of the radiative zone is considered in In §H we discuss the circumstances under 
which radiative shocks of this type are likely to occur in astrophysical sources, focusing on 
the case of GRBs and their afterglows. Our conclusions are presented in §^ 



2. The Shock Transition 



2.1. Formulation of the Problem 



We analyze the steady-state propagation of relativistic, planar, radiative shocks in a 
medium that contains a large-scale, ordered magnetic field. As we noted in §1, radiative 
shocks can be divided into a very thin "shock transition" and a "radiative zone" of finite 
thickness. We assume, for simplicity, that all the energy dissipation occurs in the shock 
transition and all the cooling in the radiative zone. The assumption of a stationary flow 
throughout the shock is justified by the fact that, by the definition of radiative shocks, the 
radiative cooling time trad of the gas that passes the shock transition is much shorter than 
the shock dynamical time tdyn- The assumption of a plane-parallel flow is easily justified for 
the shock transition, and we now argue that it is also appropriate for the radiative zone. 

The global shape of the shock is generally not planar. For example, if the shock is 
triggered by a nearly isotropic energy release at the source, as in supernova remnants 
and certain GRB scenarios, then it may assume an approximately spherical shape. The 
instantaneous distance of the shock front from the central source, R, is then given by the 
radius. If the geometry is not exactly spherical, then one needs to distinguish between the 



^Relativistic radiative shocks have previously been considered in the context of blast-wave models 

In particular, Cohen et 



( Plandford fc McKee 1976| ; pohen, Piran 
al. (1998) and Cohen 



Sari 1998| ; [Cohen fc Piran 19991) . 



Piran (1999) obtained self-similar solutions for spherical blast waves that have 
adiabatic interiors and that are bounded by fully (or partially) radiative ultrarelativistic (or Newtonian) 
shocks. These studies incorporated the efect of radiative losses in a parametric way into the overall shock 
jump conditions. They did not, however, include specific cooling mechanisms and did not consider the 
detailed structure of the shock wave. 
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distance R and the local radius of curvature Rc- Typically, Rc ^ R and both are much larger 
than the width of the radiative zone. For example, in a relativistic blast wave of Lorentz 
factor 7^1, most of the energy and swept-up matter are concentrated within a thin shell 
of width A ~ in the observer frame, or A' = 7A ~ R/'j in the shock frame (e.g., 

Blandford fc McKee 1976|) . The width of the radiative zone is ~ A(tradAdyn) -C A -C i? 



in the observer frame, or ~ A'{t[^^/t'^y^) = -Rc(^rad/^dyn) R'c the shock frame (where 
R'^ = R/j is the local radius of curvature in the shock frame). In any reasonable scenario, 
R'^ ^ A' , and even for i?^ ~ A', the width of the radiative zone would still be smaller than 
R'^ by a factor ~ trad/^dyn ^ 1- Given also that Rc is not expected to change considerably 
on time scales shorter than the dynamical time, it is thus well justified to assume that the 
shock is planar when calculating the properties of the radiative zone. 

We distinguish between the shock rest frame, in which we measure the time t, the 
distance behind the shock front x, and the fluid speed v = jSc, and the fluid rest frame 
(the frame comoving with the local mean speed of the flow), in which we measure the 
thermodynamic quantities (the particle number density n, the rest-mass density p, the 
internal energy density e, the thermal pressure p, and the enthalpy density w = pc^ + e + p) 
as well as the magnetic field amplitude B. We simplify the treatment by assuming that the 
field is purely transverse, B = B±, i.e., it only has a component that is perpendicular to 
the shock velocity and parallel to the shock front. This is justified by the fact that, upon 
traversing the shock, B_i is amplified by a factor equal to the downstream-to-upstream 
density ratio (which can be quite large in a relativistic shock) whereas the parallel 
component remains unchanged, so that the postshock magnetic field is typically almost 
completely transverse. 



2.2. Shock Jump Conditions 

We denote quantities in the unshocked upstream region (x < 0) by the subscript 1, 
whereas quantities immediately behind the shock transition are denoted by the subscript 2. 
Quantities without numerical subscripts refer to the radiative zone. The shock is assumed 
to propagate in an electrically neutral plasma consisting of protons, electrons, and possibly 
also positrons. We define the parameter 



where Up, n^- and ne+ are the number densities of protons, electrons, and positrons, 
respectively. We assume that the number of protons is unchanged by the passage through 
the shock transition, but we allow the number of pairs to possibly increase in the shock 
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transition (although not during the subsequent transit through the radiative zone). Thus 
we have 

UlUpi = U2np2 = UHp , (jyUiHel = ^2^62 = UUf, , (2) 

where Ue = n^- + nf.+ is the combined number density of electrons and positrons, u = 7/3 
is the proper speed (in units of the speed of light c), and the parameter accounts for the 
possible production of pairs in the shock. The constraint (f) > 1 implies that 

0(1 + Xi)-(1-Xi) . 
0(1 + Xi) + (1-Xi) 

From equations (P and @ we infer 

np2 = — — . 4 

\1 + Xlj U2 

With this parameterization it is straightforward to take the following limits: (i) no pair 
production at the shock corresponds to (f) = 1; (ii) no protons corresponds to Xi = X2 = 1; 
(iii) no pairs upstream corresponds to xi = 0. 

The pressure and energy density of the magnetic field are given hj ps = cb = B"^ /^-k, 
and hence the total pressure and enthalpy density are ptot = p+B'^ /8n and Wtot = w + B^/An, 
respectively. In this paper we focus on the case in which the fluid is "cold" upstream and 
"hot" downstream. The former assumption means that the thermal pressure and internal 
energy density of the upstream fluid can be neglected (ei, pi <^ PiC^), so that 

2 2 / nip 1 — Xi \ 

wi ^ pic = neimeC — — hi, (5) 

Vmel + Xi / 

where rrip and rrie are the proton and electron masses, respectively. The latter assumption 
means that the shock Lorentz factor is sufficiently large for the random motions of the 
particles of the shocked fluid to be highly relativistic, so that, everywhere behind the shock 
transition, 

Pp = ep/3, Pe = ee/3, w = 4p = 4e/3, (6) 

where pe = Pe+ +Pe-, P = Pp+ Pe, Ce = ee+ + Cg- , and e = ep + Ce- 

In the absence of a detailed shock model, it is necessary to specify how the postshock 
internal energy is divided between the proton and electron/positron components. We do 
this by parameterizing the ratio of the pressure of the electrons and positrons to that of the 
total thermal pressure just behind the shock transition, 

V = — • (7) 

P2 
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In view of the assumed postshock equation of state (eq. P), rj is just the value of the 
parameter eg (see §1) immediately behind the shock transition (77 = ee2). 

Under the assumption of ideal MHD, the magnetic flux is frozen into the fluid and 
remains constant in the shock frame.^ For a purely transverse field, this implies 

uiBi = U2B2 = uB . (8) 

The momentum and energy conservation relations across the shock transition take the form 

u\wtotl + Ptotl = tt2^f^tot2 + Ptot2 , (9) 
7l^*l^i^totl = 72M2Wtot2 • (10) 

Since our main interest is in relativistic shock waves, we henceforth set /5i ~ 1 and 
Ml ~ 7i ^ 1 and neglect the second term on the left-hand side of equation (P). Substitution 
of equations (|) and (|^) into equations (^ and (|10D then yields, after some algebra, an 
equation for /32, 



and an equation for p2, 



where 



,2„,2 ^ J2 



P"^ ~ — ^ — 2 — ' y^"^ 



Bl _ Aul + l 

iTTp2 2U2U^l — 1 



^2^7-^ = .2-2 . (13) 



is the ratio of the magnetic pressure to the thermal pressure just behind the shock transition, 
and u\ = B'^/inw is the square of the Alfven proper speed (so m^^^ = Bf/inpic^). After 
the gas enters the radiative zone it cools and gets compressed, resulting in an increase in 
the magnetic-to-thermal pressure ratio. In fact, as we show in § p.3| , 6{x) is a monotonically 
increasing function of x. 

The solution of equation ( [TT| ) is 



= ^'.2 , n • (14) 



In the limit uai ^ 1 the magnetic effects near the shock transition are weak, and one 
recovers from equations (|T^ and (|T2|) the expressions for a nonmagnetic, ultrarelativistic 



^We check on the vahdity of this assumption at the end of § 



i3.2 
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(/?! 1) shock, /?2 ~ 1/3 and p2 ~ (2/3)mJpiC^. As uai is increased, becomes larger 
while p2 becomes smaller. The dependence of (32, ^2, and 62 on uai is shown in Figure |1]. 
Note that, for a given shock speed, Ua_i has an upper bound, namely Ui. This follows from 
the fact that, in order for a shock wave to exist, the relativistic Alfven Mach number, given 
by Mai = ui/uai (e.g., [Konigl 1980|), must exceed 1. 



3. The Radiative Zone 
3.1. Radiation 

We take the radiation mechanism behind the shock to be synchrotron emission and 
inverse-Compton scattering (IC hereafter) of seed photons by the relativistic electron- 
positron component. The seed photons may consist of internally emitted synchrotron 
photons, giving rise to synchrotron self-Compton radiation (SSC), or of an externally 
produced population of soft photons, in which case the IC emission is referred to as external- 
radiation Compton (ERC). Although we treat the postshock protons as a relativistic fluid, 
we neglect their (much weaker) nonthermal emission.'' Furthermore, we neglect any possible 
transfer of energy between the proton and electron/positron components in the radiative 
zone. Under these assumptions, the protons evolve adiabatically, and the fraction of the 
energy dissipated in the shock transition that can be radiated away is determined uniquely 
by the parameter rj defined in equation (^. 

The radiating particles are assumed to have a monoenergetic energy distribution that 
can be characterized by a random (or "thermal") Lorentz factor 7e (the same for electrons 
and positrons) that is a function of the distance x behind the shock. This approximation to 
the local energy distribution is appropriate for electrons and positrons that have undergone 
significant cooling. We note in this connection that observations and particle acceleration 
models of relativistic shock sources suggest that the energy distribution immediately behind 
the shock transition has the form of a power law. In this case the radiative zone is expected 
to consist of a thin layer just behind the shock transition where most of the electrons 
and positrons have not yet cooled significantly (and are thus not well approximated by 
a monoenergetic distribution), followed by a much wider layer where substantial cooling 
has already occurred and where the monoenergetic representation is adequate (see |Granot 



I ^p ^f both the proton and electron components are effectively monoenergetic, then the ratio of their 
nonthermal (synchrotron or Thomson-scattering) cooling times behind the shock is trad,p/^rad,e — [v/i^ ^ 
" X2)/(l + X2)]('7tp/TOe)'', which is typically » 1. For a discussion of possible observable consequences 
in GRB sources of nonthermal emission by protons see, e.g., Zhang & Meszaros (2001) and references therein. 
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Piran, & Sari 2000). Although the region where the monoenergetic approximation does 
not apply is relevant to the spectral modeling of the source, it can be safely neglected in a 
calculation of the cooling-zone dynamics. 

The synchrotron emissivity (energy per unit volume per unit time measured in the rest 
frame of the fluid) from a population of relativistic electrons and positrons of total number 
density rtf. that have an isotropic velocity distribution is given by Asj^„ = aTcrie'j'^B'^ /6tt 



e.g., [Rybicki fc Lightman 1979D, where a-r is the Thomson cross section. To obtain the 



total emissivity, one must add the IC component. The IC power of a single electron is 
(4/3)(TTcne7geph, where Cph is the energy density of photons, assuming an isotropic photon 
distribution. This is also the average IC power of an electron for an isotropic distribution of 
electron velocities, even when the photon distribution is not isotropic. Therefore, the total 
emissivity is 

A = ^arcn^^j (e^ + e^^ + e^') , (15) 

where e'^^ and ep^* are, respectively, the energy densities of the synchrotron and the external 
photons in the local frame. 

We assume that the Thomson optical depth is much smaller than 1, a situation that 
often applies in astrophysical relativistic-jet sources. We further assume that the emitting 
region is optically thin to synchrotron self-absorption near the characteristic synchrotron 
frequency, which is typical in GRBs and their afterglows. Under these conditions, most 
of the radiated energy can freely escape the system.^ The average fractional energy gain, 
in a single electron scattering, of seed photons that propagate through the radiative zone 
is measured by the Compton y parameter. In the Thomson limit {"^ehv -C meC^, where u 
is the incident photon frequency in the local frame, i.e. the fluid rest frame), the average 
fractional energy gain of a seed photon moving at an angle 6 = cos~^ /i (measured in the 
shock frame) with respect to the x direction is give by 



y{f^) = ^J d'^Tll = / dx-ilne-i{l - I3^i) , (16) 



where (Itt = rieCTTds is the differential optical depth for Compton scattering and 
ds = 7(1 — l3fi)dx/fi is the differential of the distance traveled by the photon in the local 
frame. One can obtain y by averaging the expression (0) over the angular distribution 
of the seed photons in the shock frame. In the case of an infinite planar shock, the term 



®If these conditions are not satisfied, then the effective value of trad is increased by opacity effects and 
the shock may not even be radiative. In any case, if the self-absorption frequency exceeds the characteristic 
frequency, the shock will become spectrally indistinguishable from a nonradiative shock (e.g., Granot et al. 
2000). 
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l/ji that appears in equation (16) will cause the integral to diverge logarithmically for 
any angular distribution that does not vanish at /i = 0. However, under more realistic 
circumstances, one may expect the derived value of y to remain finite. We approximate the 
integral in equation (|TB|) by setting yU = 1 and then multiplying the value of the integrand 
just behind the shock transition by the characteristic size Xc of the radiative zone (see eq. 
||20| below). This leads to 

4 

y = gC^T7e2^e272(l " P2)klXc , (17) 

where we have introduced a numerical factor ki (expected to be of order unity) that 
incorporates the various uncertainties involved in our approximation. For y <^ 1 only single 
scatterings are important, whereas for y '^^1 multiple scatterings can significantly affect the 
energy and spectrum of the scattered photons. However, a typical seed photon with initial 
energy hu is upscattered on the average to an energy ~ 7e^z^ in a single scattering. Usually, 
the typical electron Lorentz factor 7e is sufficiently high for the Thomson limit to no longer 
apply even in the second scattering {'y^hu > rrieC^). This causes the energy gain in each 
successive scattering to be significantly reduced on account of the Klein-Nishina decrease 
in the cross section and of electron recoil. We therefore neglect multiple scatterings in our 
analysis. 

The energy density e^J^ of the synchrotron photons may be calculated in a manner 
analogous to that of the Compton y parameter and is subject to similar uncertainties. We 
treat it in the same manner as above by introducing a new parameter of order unity, 

4h° = ^syn— = ^(^T^e27e272(l " P2)k2Xc = ^y ( ^ ) , (18) 
CO Ki \07l J 

where x'^ = Xc7(l — (3) is the length of the trajectory of a photon in the x direction {6 = 0) 
within the radiative zone, as measured in the comoving frame. We also consider an external 
photon distribution that is isotropic and has an energy density ep^\ in the frame of the 
preshock fluid. In the local rest frame of the shocked fluid, the energy density of the 
external photons is given by 

-lih\i-pr-^lihiii-p2? (19) 




e.g., [Rybicki fc Lightman 1979| ), where 7^ = 7i7(l — l3i(3) ~ 7i7(l — /3) is the relative 



Lorentz factor between the shocked and unshocked fluids. 

The cooling length Xc is the distance behind the shock transition over which the 
electrons/positrons lose the bulk of their internal energy; it is thus the characteristic size 
of the radiative zone. Since a signiflcant fraction of the internal energy must be radiated 
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away before the hydrodynamic quantities start to deviate significantly from their values 
just behind the shock transition, a substantial portion of the radiation would be emitted 
under conditions that prevail at the upstream end of the radiative zone. One can therefore 
evaluate Xc by using the parameter values at X2, 

^ U2Cee2 _ 27Tmlc'^U2ne2 _ IGTT^m^c^nei (j)S2 , . 

A2 ~ {l + b + ky)aTBlp,2~ {l + b + ky)aTBful V ' 



where k = k2/ki and b = ^"^2/^32,^ and where we have used equations (|T5D and (0). From 
equations (pOD and ([T7|) we obtain an equation for y: 



{kyf + (1 + b)ky = k2U2l2{l -P2)— = k2U2^2{l -P2)^ = a, (21) 

eB2 02 



which yields 




(1+6)= 





1, 6^ < a 


a 


a, 6 < 1 


. a/b 


1, a < 6^ 



(22) 



In the first two limits the external photons are not important, and we obtain the familiar 
expressions for emission dominated by SSC and synchrotron, respectively ( ^ari, Narayan, fc 



Piran 19961) . The physical meaning of each of these limits can be understood by comparing 
the relative magnitudes of the different components of the energy density, keeping in mind 
that a ~ ee2/eB2, b = Cp^V^ss, and €^^2 ~ y^B2- In the first limit y ^ a^/^ ^ 1 and 6^ <^ a, 
implying that 6^2, Cp^g ^ ^ph2^ which according to equation (|T5D means that most of the 
radiated power goes into SSC. One also obtains ee2 ~ y^^Jh2, which is consistent with this 
interpretation, since, in a steady-state radiative shock, the incoming flux of internal energy 
{P2cee2 ~ cee2) cquals the outgoing radiation flux, which in this limit is dominated by SSC 
and is therefore approximately equal to cye^p^2- the second limit ?/ ~ a <^ 1 and 6 <^ 1, 
implying that e'^2^ ^ph2 ^B2i which means that synchrotron radiation is dominant over 
IC. In this case we also have ee2 ~ ^^ph2i which implies that all the internal energy goes into 
synchrotron radiation. In the third limit y a/b and l,a -C 6^, implying Cp^g, ^32^ ^'^21 
which means that most of the emitted power goes into IC-scattered external photons 
(ERC). This conclusion is consistent with the fact that ee2 ~ y^'^^2 this case, which 
means that all the internal energy is going into ERC. Note that the value of y ceases to be 



^Note that the parameter b is also approximately equal to Cp^^/esi, as both energy densities increase by 
a factor ~ 7^ across the shock transition (see eqs. |^ and ||l9|). 
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pertinent in this limit as it no longer gives the ratio of the IC and synchrotron emissivities 
(in particular, y could be <C 1 even as IC remains dominant). 

The photon energy densities Cp^* and e"^-^ are approximately constant throughout the 
radiative zone, whereas, as we demonstrate in § p.3| , the magnetic energy density cb cc B"^ 
may be considerably amplified if r/ ^ 1. Thus, even if synchrotron emission is not dominant 
immediately behind the shock transition (as in limits 1 and 3 above), it may become so 
further downstream. However, by the time this happens, only a small fraction of the initial 
internal energy would typically be left in the electron/positron component. Therefore, the 
total emission would still be dominated by IC (either SSC, as in the first limit, or ERC, 
as in the third limit). Nevertheless, given that the characteristic IC frequency is usually 
much higher than that of the synchrotron emission (by a factor ~ 7^ in the case of SSC), 
it is likely even in this case that only the synchrotron component is detected when the 
observation frequency is not too high. 



3.2. Hydrodynamics 

The steady-state momentum and energy conservation equations in the shock frame can 
be written as 

d^{u'^Wtot + Ptot) = -uA/c , (23) 
d^i^uwtot) = -lA/c. (24) 

The terms on the right-hand sides of these equations represent the contribution of 
the radiation field to the momentum and energy fiuxes, respectively. Their forms are 
appropriate when the emitted or scattered radiation exerts no net force on the fiuid; in 
particular, they apply when the radiation field is isotropic in the fiuid frame.[^ If the 
magnetic field is ordered, as we have assumed, then the synchrotron emission of any given 
particle is not isotropic in the fiuid rest frame; however, so long as the electron/positron 
rest-frame velocity distribution is isotropic, the total synchrotron radiation will carry no 
net momentum fiux in this frame and will, therefore, exert no net force on the fiuid. The 
resulting synchrotron photon energy density (eq. [0), when scattered by the isotropic 
electron/positron distribution, gives rise to an SSC emissivity that is isotropic in the 
fiuid frame. The same applies to an ERC component that arises from the scattering (by 
electrons/positrons that are isotropic in the local rest frame) of an isotropic (in the rest 



^°Note that the corresponding terms in equations (32) and (33) of Cohen et al. (1998) evidently need to 
be corrected. 
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frame of the central source) external radiation field (where now the photon energy density 
that appears in the isotropic-emission expression, in the local rest frame, is given by eq. 
ig; see |Dermer k SchlikeiserT^^ .n 



From equations (|23D and (^) one obtains an equation for the entropy: 

WtotdxU + ud^iwtot - Ptot) = -A/ c . (25) 
The magnetic field evolves according to equation (|) and the protons evolve adiabatically, 

-4/3 

(26) 

so the terms that involve them in equation ( ^5]) cancel out. One is left with 

ApedxU + 3udxPe = -A/c . (27) 

Eliminating the radiative terms from equations (|2^) and (p4|) results in 

dp -Bf^l 'iup 









(-) 







du~ 4:TTU^{1 + U'^) (1 + ^2) 

or, in terms of the dimensionless pressure p = p/p2, 

dp 262ulu~'^ — Aup 
~ (1 + m2) 

Equation (pQ]) can be integrated to give 



;i + ulf + 52 [1 - {u2/uf + 2ul Hu/U2)] 



77(1 + u^Y V \ u J 



(28) 



(29) 



P{U) = ^^^"^^ ' n 2<2 ^""^"^^"^"^^J . (30) 



Using equations (^0|), (p6l), and (0), we deduce 

~ / ^ P-{^-V)Pp .oiX 
Pe M = = 31 

(1 + n2)2 + ^2 [1 - (m2/m)2 + 2ul\n{u/u2)] 1 - ^^2^/^ 



^^We ignore in our discussion the effect of the radiation drag force (corresponding to a nonvanishing 
momentum flux of intercepted photons in the fluid rest frame) that is exerted by an external radiation field 
on the bulk motion of the fluid. This is consistent with our neglect of any variations in the gross properties 
of the shock that occur on time scales ^ tdyn- 
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where Pp = Pp/pp2 and Pe = Pe/Pe2- We already have exphcit expressions for 
Pe, Pp, B, Up, Tie, and 7e = 3pe/ne'meC^ as functions of u. To obtain their values as 
a function of the distance behind the shock, we need to know u{x). Introducing the 
dimensionless variable C, = x/xc and using equations (|T^, (|T3|), (PU]), (P?]), and we get 



du 
d^ 



37]uluil + u^)piu)l [(1 + Y)/{1 + Y2)] 



4(1 — 2^2) {1 — 7])u2^^u'^/^ + riu'^p{u) 



Q62U2 



where the function 



Y(u) 



(b + ky) (-) 

\U2J 



(32) 



(33) 



represents the ratio of the local IC (including SSC and ERC) and synchrotron emissivities, 
and Y2 = {b + ky) is the value of Y just behind the shock transition. Note that the only 
dependence of the function Y{u) on the Compton y parameter (see eq. [|17|) is through the 
combination ky, which, in turn, is given (in terms of the parameters a and b) by equation 



(p^). As both U2 and 62 are functions of uai (see eqs. |T3| and |T^), a is determined by the 
parameters /c2, mai, and r]. For simplicity, we set ^2 = 1 in the rest of this paper. Therefore, 
Y{u) and the right-hand side of equation ( P^ are specified by the parameters rj, uai, and 
b, and m(^) can be obtained by solving equation (^2|) numerically, using the boundary 
condition = 0) = U2- 

Before proceeding to present the solutions of the above set of equations, we digress 
briefly to check on the validity of the ideal-MHD approximation that underlies our analysis. 
Under the assumption of ideal MHD, the electric field vanishes in the rest frame of the 
fiuid. The magnetic field in the shock frame is then given by = uiBiCy/ (3{x), where we 
have chosen the y axis to point along the direction of the magnetic field. This implies a 
current density 

UiBic dp ^ 



J 



4tt 



V X (75) 



(34) 



47r/?2 dx ■ 

Since in our case j is perpendicular to the direction of the velocity, it has the same 
magnitude in the shock frame. In order for the ideal-MHD approximation to be valid, the 
current density in the fiuid rest frame cannot exceed j^ax = nec = n2U2ec/u{x), where 
n = He + Up is the total number density of charged particles (e.g., [Melatos fc Melrose 
1996|) . Furthermore, so long as j <^ jmax, the effect of anomalous resistivity associated 



with current-driven plasma instabilities is expected to be small (e.g., [Spruit, Daigne, fc 
Drenkhahn 2001| ). The condition j < jmax translates into 



7 
/3 



d/3 
dx 



< 



4'7ren2 

'b7~ 



(35) 
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In the radiative zone / I3)\d(3 / dx\ ^ 72/a;c ~ so equation (|35|) may be written as 

'i'Ken2Xc 64:7i^mlc^e nineiU2 f n2U2\ 



Note that this inequahty can be rewritten in terms of the postshock electron Larmor radius 
rL,2 as rL,2 < 3xc/262. This shows that, if ideal MHD holds in the radiative zone, then 
the condition for the fluid approximation to be valid there (essentially rL,2 <^ Xc) will also 
be satisfied. We now check whether the inequality ( |5BD holds under the most unfavorable 
choice of parameters. The latter is obtained by taking 77 = = 1 and ^2 1, which implies 
n2U2 = riiUi, uai <^ 1, M2 ~ 1/v^, and 82 ~ 3Bl/2npic^, resulting in 

2aT{l + Y2)p,Bful ■ ^'"^ 

The minimal value of ngi/pi is l/m^, which is obtained for xi = 0. This leads to 

3n^mye _ 9.5 x lOSpi^p 2.2 x 10%i,o 

^ 2aTm,(l + Y2) Bfuf " (1 + Y2)Bl,ul (1 + ^2)51,0 ' ^ ^ 

where here and below numerical subscripts (such as the subscript in ripi 0) refer to powers 
of 10 in c.g.s units. 

As an illustration we apply this condition to GRBs, considering both the prompt 
high-energy emission (resulting from internal shocks within the flow) and the lower-energy 
afterglow (arising from the interaction of the flow with the surrounding medium). We 
begin with internal shocks. We use the conventional parameterization = 8TreB2G2 (where 
we identify €32 = ^2/3), and write (e.g., Piran 1999) 62 = (712 — l)(47i2 + 3)npimpC^, 
Upi = E / iTrnipR'^'y'^c^tduT, and R = 27^ctdurC! where E is the total energy deposited in the 
flow, 7 is the typical Lorentz factor of the flow with respect to the central object (and the 
observer), 712 is the relative Lorentz factor between the shocked and the unshocked fluids, 
tdur is the duration of the burst, and ( = St/tdur is a measure of the variability of the burst, 
with St being the smallest time scale for signiflcant variation during the burst. Thus we 
obtain 

^ ^ 1.3 X 10^t^/^,,,7lC-2 ^ 3.6 X 10^t:^/^,,7fC-2 

(1 + ^2)^2^5^ [(7i2-l)(47i2 + 3)f/' ^ il + Y2)e'£Ell' ' 

where we used 712 = 2 on the right-hand side (consistent with the expectation that 712 is 
of the order of 1). Turning now to the afterglow stage of GRBs, we evaluate the inequality 
(W^) by using 62 = ij'^npimpC^ in the parametrized expression for B2. This yields 

1< 1.2 X 10^ (1 + Y2)-'e],fn;,fj2' , (40) 
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where we scaled 7 and n^i by their rough upper hmits for typical afterglows. (The Lorentz 
factor 7 decreases during the afterglow evolution, as does Upi if the shock propagates into 
a stellar-wind environment.) Since the factor (1 + Y2)e^^2 ™^ equations ( P^D and (^D|) never 
exceeds 1 in the absence of a strong external radiation field (being approximately equal to 
max{e^4^, eB2??^^^}; see eq. |^), we conclude that the use of the MHD approximation for 
studying the cooling layers behind radiative shocks should be well justified throughout the 
evolution of typical GRBs and their afterglows. 



3.3. Results 

At large distances behind the shock (,^ ^ 1) the electrons and positrons are left with 
only a small fraction of their internal energy, and all hydrodynamic quantities approach 
some constant asymptotic value: u ~ Wmin, B ^ -Bmax = -B2W2/Mmin- By equating the 
radiative cooling time with the flow time, one infers that 7e ~ 7e2/^ in that region (e.g., 
Granot et al. 2000). Also, Ue ~ ne2M2/wmin, and hence pe ~ M2/Mmm^ for ^ ^ 1. The fact 
that peiO > for any finite value of ^ implies that u{^) > Umin, so u becomes equal to Wmin 
only asymptotically. One can determine -Umin by setting equation ( |3TD equal to zero and 
(numerically) solving for u: this yields Mmin as a function of r] and uai- In Figures |^ and ^ 
we fix one of these variables and show Wmin, normalized by M2, as a function of the other 
variable. In the limit 77 — > there are no radiation losses, so the speed of the shocked fluid 
does not change: Mmin = M2- In the limit 77 — *■ 1 the electrons and positrons contain all the 
postshock internal energy and radiate it away, and the attendant compression eventually 
renders the magnetic pressure dominant. Since the magnetic field is amplified by a factor 
~ U2/umin in the radiative zone (see eq. 0), a smaller initial magnetic field (smaller uai) 
results in a lower value of Mmin/M2- Furthermore, the total pressure, ptot, is approximately 
constant throughout the radiative zone, and for uai <^ 1 it becomes independent of uai 
(Ptot ~ P2 ~ (2/3)'U^pic^). Since in the limit rj 1 the magnetic pressure eventually 
becomes dominant (-Bmax/^Ti" ~ Ptot), -Bmax also becomes independent of uai for uai <^ 1, 
and we obtain Umm/u2 = -B2/-Bmax oc uai, as can be seen in Figure ||. 

Figure ^ displays u{C,), normahzed by M2, for 77 = 0.9, b = 0, and several values of uai- 
For large values of uai {'^ 0.08) the asymptotic value u = Umm (shown by a dashed line) is 
approached at smaller values of ^ (traced by the asterisks) as uai decreases, whereas the 
converse occurs for smaller values of uai- This can be understood as follows. For large 
values of uai the initial magnetic field B2 is strong and Y2 is close to 1, so synchrotron 
cooling is important early on. As the internal energy is radiated, the fluid becomes denser 
and the magnetic field grows stronger, which enhances the synchrotron emission. Since 
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the cooling length is calculated by using the initial synchrotron emissivity, the value of Xc 
is overestimated in this case and hence the asymptotic value of u = U2 is approached at 
lower values of ^. This effect is stronger the larger the value of -Bmax, which for a fixed 77 
is obtained for lower values of uai- However, for sufficiently small values of uai, Y2 comes 
to exceed 1 and IC becomes the dominant cooling mechanism. For uai <^ 1 synchrotron 
radiation hardly affects the dynamics, and Xc is determined by the IC emissivity just 
behind the shock transition. Since, unlike the magnetic energy density, the photon energy 
density that determines the IC cooling rate does not increase with ^, the value of Xc better 
approximates the actual coohng length, and consequently Wmin is approached at larger 
values of ^ for smaller uai- 

Figure ^ shows normalized by M2, for uai = 0.05, 6 = 0, and several values of 
7]. In this case the asymptotic value of u is approached at smaller values of ^ for larger 
values of r/, which correspond to larger values of -Bmax and a stronger enhancement of the 
synchrotron cooling. 



It is straightforward to verify from equations (^) and (|3^ ) that Pe is a monotonically 
increasing function of u and that m is a monotonically decreasing function of x. Therefore, 
as expected, both pe and u decrease the further one gets from the shock front. 

To demonstrate the evolution of the various components of the pressure, we show in 
Figure ^ the dependence of Pe, Pp, Vb and ptot (each normalized by P2) on m for 77 = 0.9 
and Uai = 0.2. Equipartition occurs when the total thermal pressure equals the magnetic 
pressure, p{ueq) = psiueq), with Meq denoting the value of u at equipartition. However, this 
condition is not satisfied in all cases. There is a critical value of uai, which we denote by 



eq; 



for which 62(u^ 



1. We find that Ucq = 0.447. For uai > u^q the magnetic field just 
behind the shock transition is already above equipartition, B2 > -Bcq2; where B^q is given 
by -Beq = (87rp)^/^. On the other hand, B2 < -Beq2 for uai < Ueq, and only if t] is above some 
critical value rj^miuAi) will the asymptotic magnetic field be above equipartition. Since U2, 
and therefore also 62, are functions of uai, p is a function only of uai and u (see eq. [0)- 
Therefore, 77min(wAi) may be found by (numerically) solving the equation 



p[uAl,Unim{uAl)] = S2{uai) 



U2{uai) 
UminiuAl 



(41) 



where Umin('"Ai) = u^m[Vnim{uAi),UAi]- Figurcs and ^ show r^min as a function of uai and 
82, respectively. The critical value of rj is close to 1 for small values of uai (or ^2) since, in 
this case, the initial magnetic field is small, and almost all of the internal energy needs to 
be lost in order for -Bmax/-B2 (= U2/umm) to be large enough for -Bmax to reach equipartition. 



As we noted in 51, there are direct observational indications that at least in some 
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relativistic shock sources the magnetic field is not far below equipartition. Using the 
calculated radiative shock structure, we can address the question of how far from 
equipartition the magnetic field is at the location where the bulk of the observed radiation 
is emitted. A fiuid element starts radiating just behind the shock transition, where B = B2 
and 6 = 62, and as it fiows downstream the magnetic field increases, and so does the 
magnetic-to-thermal pressure ratio 6. Thus, even if initially the magnetic field were much 
below equipartition (^2 <^ 1), it would be amplified by a large factor (5 ^ ^2) in the 
radiative zone and might eventually even exceed equipartition (5 > 1). It is therefore 
conceivable that a considerable fraction of the radiation could be emitted from a region 
where the magnetic field is much larger than the immediate-postshock value. We now 
examine this possibility in a quantitative manner. 

Since the internal energy of the electrons and positrons is characterized by their 
thermal Lorentz factor 7e oc Pe/^e oc upe-, their fractional remaining energy is given by 
upeiu) /u2Pe2 = upe{u)/u2. Figure P shows this quantity as a function of 82 at two different 
instants: (i) at equipartition, when 5 = 1, and (ii) when 5 = 0.01. For rj close to 1 a 
significant fraction of the internal energy of the electrons and positrons is still left when 
the field reaches equipartition even if B2 -C -Beq2- In this case a significant portion of the 
radiation originates in a region where the magnetic field is close to equipartition. We find 
that, so long as the energy dissipated in the shock transition is deposited mainly in the 
component (i.e., 77 1), then, even if 62 is as small as ~ 10~^ — 10~^, a considerable fraction 
of the radiation is still emitted at 5 > 10^^. Note that, for 77 1 and ^2 ^ 1, a ~ 1/^2 ^ 1 
(see eq. [^), so that, in the absence of a strong external photon field (6 ^ a^^^ ~ ^2), 
y ^ a^/^ ~ ^2 ^ 1 (see eq. [^). This implies that SSC dominates the synchrotron 

1 /2 

emission, with the synchrotron component contributing only a fraction 1/(1 -(-|/) ~ ^2 1 
of the total radiation. However, as we noted at the end of §^[1] (see also §|), the synchrotron 
emission may dominate the SSC contribution in the observed frequency range, in which 
case a significant fraction of the observed radiation could, in fact, originate from a region 
where S ':$> 62- 



In Figure |10| we show the fraction of the total synchrotron emission that is emitted 
above a given value of B/Beq, as well as the fraction of the total emission contributed by 
the synchrotron process above that value of B/Beq, as a function of the magnetic field 
strength behind the shock. Three sets of parameters are shown, labeled A through C, 
which correspond to the three limits (from top to bottom, respectively) of equation (0). In 
case A, SSC is dominant near the shock front, whereas further downstream the magnetic 
field is strongly amplified by compression, and synchrotron emission becomes dominant. 
Altogether, synchrotron emission constitutes 45% of the total emission, whereas near the 
maximal field amplitude {B ^ Bj^so^ the synchrotron process contributes 94% of the total 
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emission. A large fraction of the synchrotron radiation comes from a region where the 
magnetic field is close to equipartition even though the magnetic field just behind the 
shock transition is a factor of 10 below the equipartition value (^2 ~ 10~^). In case B, 
synchrotron emission is dominant and contributes 84% of the integrated total emission. 
The compression behind the shock is small and the magnetic field is hardly amplified, so 
even near -Bmax the synchrotron fraction of the total emission (86%) is only slightly higher 
than the overall fraction. In case C, the values of rj and uai are the same as in case A, and 
therefore the initial and final values of the magnetic field are the same. However, in case C 
the emission is dominated by IC scattering of external photons {b — 100), which are absent 
in case A. Synchrotron emission now constitutes only 13% of the integrated total emission, 
but near -Bmax it contributes 59% of the total. Most of the synchrotron emission comes from 
a region where the magnetic field is close to equipartition — a larger fraction than in case 
A. This is a reflection of the fact that the ratio of synchrotron to total emission near i?max 
relative to that near B2 is larger in case C because of the added cooling by IC scattering of 
the external photons. 



4. Applications 

To assess the relevance of cooling-induced magnetic field amplification to relativistic 
shock sources, one needs to address the following questions: What are the likely values 
of the normalized energy densities of the electrons/positrons and the magnetic field (cg 
and es, respectively; see §1) in the emission region? What are the expected values of 
and €3 immediately behind the shock transition (which can be written, in terms of the 
shock parameters defined in §2, as ee2 = rj and €32 — ^2/^)'^ Is the shock radiative? If 
the shock is radiative but €32 is lower than the observationally estimated value of e^, can 
the compressional amplification of the magnetic field in the radiative zone make up the 
difference? 

A strong relativistic shock (i^i ^ 1, uai <S 1) with a transverse magnetic field satisfies 
62 ~ 2ulwi and S| 8ulBf (see §2), implying that the ratio of the magnetic energy 
density to the total matter energy density is approximately constant across the shock and 
hence that €32 can be expressed in terms of the enthalpy density Wi (= pic^ for a "cold" 
upstream medium) and the magnetic field amplitude Bi ahead of the shock: 



As a concrete example of how the above considerations can be applied, we again focus 
on GRBs and their afterglows. However, at the end of this section we comment briefiy 
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on AGNs. Starting with the 7-ray bursts themselves, we have aheady observed (see §1) 
that efficiency considerations imply that eg (in particular) and cannot be much smaller 
than 1 in the emission region. The value of may be constrained by observations through 
the requirement that the cooling time does not exceed the variability time scale 5t. This 
translates into the condition that the Lorentz factor of an electron that cools on the time 5t, 
67rmec/aTB'^'j6t{l + 1^2)5 must remain smaller than the characteristic synchrotron Lorentz 
factor, (27rmeCz/ob/7e-B)^/^, that corresponds to the observed frequency z/ob-Q Following the 
same steps as in the derivation of equation (^) and again setting 712 = 2, we deduce 

> 8.3 X 10-(1 + F2)-/^i^5-.^72"/^t^/liC!? (j^) (43) 



(cf. [Sari fc Piran IQQTal) . When the energy density of external photons is negligible, one has 



Y2 ^ y (see eq. with b = 0), so that, ii y ^ 1 (in which case y ~ {ee2/^B2y^'^', see eqs. 
2TI and [^), the constraint (^3|) on would be relaxed. However, y ^ I means that IC is 



dominant, with the observed synchrotron radiation being only a fraction ~ l/y of the total 
radiated energy. The requirement of a reasonable radiative efficiency would then constrain 
es to be not much smaller than eg (e.g., Piran 1999). The latter restriction would not arise 
if the observed 7-ray emission represented the SSC rather than the synchrotron component, 
but the synchrotron model appears to account quite well for the observed spectra (e.g., 
P^loyd fc Petrosian 2000|) , whereas an SSC interpretation could entail a prohibitively large 



energy input at the source .0 

As we have already noted, the prompt 7-ray emission in GRBs is believed to be the 
result of internal shocks within the flow. The relevant preshock field is then the relic 
magnetic field advected from the origin, and, to the extent that flux freezing is maintained. 



^^To simplify the discussion, we omit cosmological correction factors from the expressions presented in 
this paper. 

^•^For example, in the model of Panaitcscu and Meszaros (2000), the observed 7-ray spectrum is identified 
with a singly scattered SSC component. The characteristic thermal Lorentz factor of the emitting electrons 
in this picture is ~ 30, which implies that the peak synchrotron frequency is not very far from the optical 
regime. There should thus be an optical pulse that coincides with the prompt 7-ray emission, and the fact 
that such a component has not been detected (e.g., Akerlof et al. 200C ; Williams et al. 2000| ) indicates that 



the Compton y parameter is very large (Panaitescu and Meszaros 2000 estimate y ^ 5000). Now, for the 
postulated thermal Lorentz factor, a typical seed photon undergoes two scatterings before the Klein-Nishina 
limit is reached. Since the emissivity ratios of the synchrotron, singly scattered SSC, and doubly scattered 
SSC components are 1 : y : y^, it follows that the observed 7-ray pulse may constitute only a fraction 
~ l/y ^ 1 of the total burst emission (and hence that the energy requirement on the source is a factor of y 
larger) . 
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eB2 oc Bf/pi keeps the value that was imprinted on the ejecta at the source.[^ Spruit et 
ah (2001) pointed out that, under these circumstances, values of eB2 as high as ~ 0.1 — 1 
may be expected in the commonly invoked GRB models that envision large-scale, ordered 
magnetic fields actively driving the relativistic outfiow by tapping the rotational energy of 
the central source. In particular, €32 could well be ~ 1 in coUimated (jet-hke) outflows of 
this type. Spruit et al. further estimated that, if the fleld is "passive" (in that it does not 
drive the flow but is merely advected outward from the source), then 6^2 would lie in the 
range ~ 10"^ - 10~l In the latter case, compressional magnetic fleld ampliflcation in the 
radiative zone could increase the fraction of radiation going into the synchrotron component 
(see Figs. ^ and |10]). For this mechanism to operate, it is necessary that the radiative 
cooling time be much shorter than the dynamical time and that eg be ~ 1. The flrst 
requirement should in general be satisfled in internal shocks that produce GRBs (e.g., Piran 
1999). The conditions under which eg could be close to 1 are somewhat uncertain, but there 
are several conceivable mechanisms for depositing a signiflcant fraction of the postshock 
internal energy in the electron component even if the outflow consists predominantly of an 
electron-proton plasma (e.g., Pykov fc Meszaros 1996|) . It should be even easier to attain 
high values of eg if the plasma is dominated by pairs, a situation that in fact arises 
naturally in some GRB scenarios that involve magnetically driven outflows (e.g., Usov 1994; 
IMeszaros fc Rees 1997|) . 

The long-lived afterglow emission of GRBs is most likely associated with the forward 
shock that propagates into the ambient medium. We have already noted in §1 that 
comparatively high values of eg (~ 0.1 — 1) and (~ 0.01 — 0.1) have been directly inferred 
from observations of the GRB 970508 afterglow. These parameters have so far been less well 
constrained in other sources, but in a recent modeling of four other afterglows ( |Panaitescu 
fc Kumar 2001| ) values in the range 10~^ ^ Ce ^ 10~^ and 10~^ ^ es ^ 10~^ have been 
inferred. It is also worth noting in this connection the suggestion by Galama et al. (1999) 
that comparatively low values of may characterize afterglow sources that exhibit weak 
radio emission (attributed by them to low values of the characteristic frequency and peak 



^ Note in this connection that equation ( [42[ ) shows quite generaUy that €32 wiU remain constant with 
time if the medium into which the shock propagates is either a supersonic outflow with a predominantly 

1 /2 

transverse magnetic field that is frozen into the matter (in which case both Bi and Pi scale inversely with 
distance from the center) or else is uniform (in which case the values of Bi and pi are fixed). The former case 
applies to internal shocks and to the reverse shock (which propagate in the ejecta) as well as to a forward 
shock that propagates in a magnetized stellar- wind environment. The latter case corresponds to a forward 
shock in a uniform ISM. This shows that the predictions of conventional emission models of GRBs and their 
afterglows, which simply postulate a constant postshock value of e^, will by and large continue to apply 
when this parameter is calculated self-consistently. 
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flux of tlie synclirotron spectrum, botli of wliicli scale as ). The expected values of ee2 
and eB2 in GRB afterglows depend on the GRB environment (which directly affects such 
factors as the pair content, the electron/positron acceleration efficiency, and the degree of 
magnetization in the shocked gas). The uniform-density model fittings of the GRB 970508 
afterglow ( [Wijers fc Galama 1999| ; |Granot, Piran, fc Sari 199UD yielded upstream densities 



in the range ~ 0.03 — 3 cm"'^, which are typical of a diffuse ISM. In the case of GRB 980329, 
a preshock density ~ 10^ cm~^, which is typical of a molecular cloud, has been deduced 
( [Lamb et al. 199"9D . In the Galactic ISM, the inferred magnitude of the Alfven speed does 
not vary strongly from diffuse to dense environments and corresponds to values of 6^2 in 
the range ~ 3 x 10^^^ — 3 x 10^^°. If these values are also typical of the host galaxies of 
GRBs, then the synchrotron emission from their associated afterglows would in most cases 
be too weak to be detected. However, all the afterglow sources observed so far correspond 
to "long" bursts (tdur ^ 2 s), for which a plausible origin is the collapse of a massive 
(e.g., Wolf-Rayet) star (e.g., Woosley 2000; Meszaros 2000). In this case, the medium into 
which the forward shock expands could correspond to the preburst wind from the central 
star (e.g., [Chevalier fc Li 2000t |R amirez-Ruiz et al. 2001[) . It has been suggested (e.g.. 



Biermann fc Cassinelli 1993[) that winds from Wolf-Rayet stars are magnetically driven and 



have spatially constant Alfven speeds corresponding to values of €32 that can be as high as 
~ 10~^ . If such winds constituted the environment into which the forward shock expands, 
then they could account for afterglows with the lowest values of inferred to date even if 
there were no additional magnetic field amplification beyond that induced by the passage 
of the ambient field through the shock transition. Radiative cooling effects could increase 
the magnitude of in the afterglow emission region above its value immediately behind 
the forward shock. The cooling time behind the forward shock is typically shorter than the 
dynamical time during the early phase of the afterglow, but subsequently this inequality is 
reversed. For a shock propagating into a uniform- density ISM, the time of transition from 
radiative to adiabatic evolution was calculated by Sari et al. (1998) and is given by 

trad^ad = 4.6 e'£eli'{l + Y,y/'Eti'ioi^'nT days (ISM) , (44) 

where 70 is the initial Lorentz factor of the ejecta. In the absence of a strong coupling 
between protons and electrons behind the shock transition, ee2 ~ 1 is required for the shock 
to be radiative. Equation (^^ indicates that, even if this condition is satisfied, (1 + ¥2)632 
cannot be much smaller than 1 if the radiative phase is to last through a significant fraction 
of the early evolution of a typical afterglow (i.e., a few days or so). Taking 6^2 ~ 10~^° as a 
representative ISM value, we infer (1 -t- ¥2)632 ~ {^€2^32)^^'^ ~ 10~^, implying that trad^ad 
is only a few seconds! In fact, in this case, even if the radiative phase lasted longer, the 
expected magnetic field amplification in the radiative zone would not be sufficient to raise 
to a detectable level (see §3.3[). 
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In the case of a stellar-wind environment, the ambient density scales with distance 
r from the source as p(r) = Ar^'^, where the constant A = M/AttVw is defined in terms 
of the wind mass-loss rate M and speed f^. We are interested in the evolutionary phase 
during which the forward shock has already decelerated significantly but is still relativistic. 
Assuming that the mass is concentrated in a thin, spherical shell of radius R and Lorentz 
factor 7, we can write 7(-R) = M/m{R), where M (= E/jqc^) is the initial mass of the 
ejecta and m{R) is the rest mass of the swept-up matter QBlandford fc McKee 1976] [Katz 



fc Piran 1997|) . In our case m{R) = AttAR and therefore 7(-R) = L/R, where L = M/AnA 



is the distance from the center where the flow becomes nonrelativistic. Using t = R/4'y'^c, 
we obtain 



167r" 



-Ht.^ . (45) 



/ \ 1/3 

The Lorentz factor of an electron that cools on a time scale t is given by 

_ QirrrieC 

whereas the typical Lorentz factor of the electrons can be approximated by 

7^ ^ 610ee27 (48) 

(e.g.. Sari et al. 1998). Assuming that the preshock field Bi = -B2/71 scales as 1/i? (e.g., 
piermann fc Cassinelli 1993|) , we find 7c(t) oc t"^^^, so the ratio of the radiative cooling time 
to the dynamical time (trad/^dyn = Ic/lm) is oc t. This shows explicitly that at early times 
^rad < ^dyn and the evolutiou is radiative. The transition from the radiative to the adiabatic 
phase occurs approximately when ^dlm increases to 1. Using B"^ /Sn = e^e and e = 47^pc^, 
we obtain 

trad^ad = 28.8eB2ee2(l + ^^2)^* days (Wind) , (49) 

where A^ = A/{5 x 10^^ g cm~^) corresponds to a stellar wind characterized by a speed of 
10^ cm s~^ and a mass loss rate of 10~^ Mq jt~^ ( [Chevalier fc Li 2000|) . As can be seen 



from a comparison of equations (|4^) and (^), the radiative phase of the forward shock will 
typically last longer in a wind environment than in a uniform- density medium, which can be 
attributed ( [Chevalier fc Li 2000|) to the higher ambient density encountered at early times 



by a shock that propagates into a wind. Using equations ([2l[) and (p2[) to approximate 
^e2^B2{^ + ^2) ^ (e_B2/2)"'^^^ and setting €32 ~ 10""^ (the expected value in a strongly 
magnetized Wolf-Rayet outflow), we infer from equation (^) that trad^ad ^ SA* hr. During 
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the radiative phase, a typical electron could radiate up to ~ 10% of its postshock random 
energy in a region with ;^ 10^^ (see Fig. ^. However, the estimated value of trad^ad is 
too small for the radiative field-amplification mechanism to apply to the observations of a 
source like GRB 970508, in which values of ^ 10~^ have been inferred over substantially 
longer time scales. Furthermore, even if the model were applicable, this mechanism could 
not account for a value of as high as 0.1 (the best-fit value for a stellar- wind model; 
Chevalier fc Li 2000|) if €32 is not larger than ~ 10~^. 



Another potential source of detectable emission in GRB sources is the reverse shock, 
which runs into the ejecta shell when the latter starts to be decelerated by the inertia 
of the swept-up ambient gas. In fact, a shock of this type is believed to be the source 
of the prompt optical emission (or "optical flash") that was observed in GRB 990123 
( [Akerlof et. al. 1999| ; Sari & Piran 1999a,b; [IVleszaros fc Rees 1999| ) as well as of a rapidly 



decaying "radio flare" detected about a day after the burst in this source (and evidently 
also present in GRB 990123; Kulkarni et al. 1999). So far, however, these observations 
have not provided strong constraints on the parameters eg and in the emission region. 
It may, however, be expected that, since the reverse shock propagates into the ejecta, it is 
characterized by values of and €32 thai are similar to those in internal shocks. 

GRB 990123 was a comparatively long 7-ray burst: such sources likely possess a 
relatively thick ejecta shell (of thickness A as measured by a stationary local observer), and 
their prompt optical emission should peak at tpeak ~ ^/c, the shock crossing time of the 
shell (which is approximately equal to the duration of the burst), by which time the reverse 
shock would be relativistic (e.g., ^ari fc Piran 1999b| ). This is a natural time scale to use in 



a comparison with the cooling time for determining whether the shock is radiative. We can 
repeat the procedure applied above to the forward shock by comparing the characteristic 
Lorentz factors 7^ and 7c, except that in this case we substitute for 7 in equation ( ^7] ) 
for 7c, by the Lorentz factor 73 of the shocked shell material, which is given in terms of 
the Lorentz factor 7sh of the un-shocked portion of the ejecta shell (both measured in the 
stationary frame) by 

73 = %l"f^VV2 (50) 



e.g., |Sari fc Piran 1995|) , where 



E 

^ ATi'^l^R\AnpimpC^ ^^^^ 

is the ratio of the proper mass densities in the shell and the ambient medium (assumed to 
be dominated by protons) and -Ra = 73 A is the radius where the reverse shock crosses the 
shell. Similarly, we replace 7 in equation (^) for 7^ by the relative Lorentz factor between 
the shocked and the unshocked shell material, 73 = 73/a/7- 
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In the case of a uniform ambient medium, the condition 7c(tpeak)/7m(^peak) < 1 imphes 
ee2eB2(l + Y2) > 0.12 E52'^"n;f/'(tpeak/50 s)-^/^(7sh/200)-^ , (52) 



where we have normahzed tpeak and 7sh by their estimated values in GRB 990123 ( |Sari| 
fc Piran 1999a| ).p| Assuming ee2 ~ 1 ^ [which imphes, in the absence of an external 



radiation field, 662652(1 + Y2) ^ (£52/2)"*^^^], we obtain a lower bound on 6^2, 

eB2 ^ 0.03 E52'/'n;f /'(tpeak/50 s)-i/2(^^^/200)-2 . (53) 

The inequality ( |55D could in principle be satisfied if the GRB outflow is driven magnetically, 
in which case the ejecta could be characterized by e^i ~ 6^2 ~ 1- However, even if €32 
is near its estimated lower bound, compressional amplification could boost the magnetic 
field to equipartition values in the radiative zone (see Figs. ^ and P). The corresponding 
condition for a stellar wind environment (assuming strong hydrogen depletion as appropriate 
for a Wolf-Rayet progenitor) is 

ee2eB2(l + Y2) > 9.6 x IQ-^EU^A^'/^tp^^JbO sY/\%i,/200)~' (54) 

( phevalier fc Li 2000| ). For ee2 ~ 1 ^ 6^2 and a negligible external radiation field, equation 
(p^) implies 

eB2 1-8 X 10-^°^52A:'(Wk/50 s)(7sh/200)-' , (55) 

which demonstrates that the reverse shock is always radiative in this case for plausible 
parameter values. It is worth noting, however, that if the reverse shock is radiative at the 
time that it crosses the shell, then it will rapidly cool and become invisible at later times 
(e.g., |Sari fc Piran 1999a|) . Since this was not the case for the optical flash in GRB 990123, 



we can conclude, on the basis of equations ( ^4]) and (|52|), that the outflow in this source 
does not expand into a wind (as was already inferred by [Ghevalier fc Li 2000|) and that the 
value of ee2^B2 behind the reverse shock is < 0.1. 

We can summarize our discussion of GRB sources as follows. The field amplification 
mechanism discussed in this paper could in principle operate in shocks that are driven 
into the ejecta: either the internal shocks that give rise to the prompt 7-ray emission. 



^^The condition 7c = 7m is probably a bit conservative. For one thing, as we have demonstrated in §3, 
the cooHng induces compression in the radiative layer, which has the effect of enhancing the radiative losses. 
Furthermore, the radiative efhciency decreases only gradually after 7c comes to exceed 7„i: in fact, the 
fraction of the internal electron/positron energy that is radiated by the synchrotron process after that time 
is given by (7m/7c)^~^, where p (which has a canonical value of 2.5) is the power-law index of the radiating 



particle distribution (e.g., Moderski, Sikora, & Bulik 2000). Thus, the right-hand sides of the inequalities 



(p4)-(pq) are probably only upper limits. 
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or the reverse shock that produces the optical flash and radio flare emission. Except for 
reverse shocks in sources with a uniform ambient medium in the case when the condition 
(0) is not satisfied, these shocks will be radiative if most of the energy dissipated in 
the shock transition is deposited into the electron/positron component. Compressional 
field amplification in the radiative zone of these shocks could enable a significant fraction 
(~ 0.1 — 0.3) of the radiated energy to be emitted from regions where the field is substantially 
stronger (by up to ~ 1 — 2 orders of magnitude in e^) than immediately behind the shock 
transition. If the preshock magnetization is comparatively weak (e^i ^ 1), the radiative 
amplification could significantly enhance the synchrotron emission and may in certain cases 
help to bring it above the detection limit. If the ejecta is already moderately strongly 
magnetized (e^i ^ 10^^), this process could increase the field in the emission region to 
equipartition values. 

We have also concluded that the compressional amplification mechanism will not 
significantly affect the afterglow emission of GRB sources, both because the forward shock 
typically undergoes a transition to adiabatic evolution early on, and because, even during 
the radiative phase, the expected amplification falls short of the required enhancement for 
standard ISM environments. We noted that magnetized stellar winds from Wolf-Rayet stars 
could provide an adequate environment in which afterglows with values of that lie near 
the low end of the observationally inferred range would be produced in adiabatic shocks 
without requiring further field amplification. However, the comparatively high 10^^) 
values of indicated in a source like GRB 970508 cannot be accounted for by the basic 
shock model and typical ISM parameters. We do not pursue this topic further in this paper, 
but we note that several suggestions that bear on the need to have high values of and 
ee in such shocks have already been made in the literature (e.g., [Medvedev fc Loeb 1999|; 



Smolsky fc Usov 2000| ; I'l'hompson fc Madau 2000| ; [Lyutikov fc Blackman 2001| ; |Konigl fc 



Granot 20011) . 



The observational data on spectral flares in AGNs are not as extensive as in GRB 
sources, and, correspondingly, their theoretical study has so far been less developed. 
However, as we noted in §1, the blazar flares have been similarly interpreted in terms 
of nonthermal emission from relativistic, and likely highly collimated, outflows from the 
vicinity of a central compact object.|3 In particular, the variable emission has been modeled 
in terms of internal shocks in a magnetized relativistic outflow, possibly formed from 



16 



One key difference between tfie fiigli-energy flares in AGNs and GRB sources is in tfieir observed duration, 
which is of the order of liours in blazars (e.g., Wehrle et al. 1998 ; Maraschi et al. 1999 ) as compared with 
seconds in GRBs. A related difference is in the inferred bulk Lorentz factor of the outflow: < 10 in blazars 



vs. > 10^ 



GRBs. 



- 28 - 



overtaking collisions of disturbances ("shells") in the flow (e.g., [Marscher fc Gear 1995 



Romanova fc Lovelace 19"97| ; [Sikora et al. 1997| ; Levinson 1998, 1999). In these models, the 
high-energy emission is attributed to an SSC or an ERG component. The quasi-isotropic 
external radiation field invoked in the ERG interpretations can arise from the nuclear 
accretion disk (e.g., permer fc Schlikeiser 199^) , or, perhaps more likely, from scattering 
by the broad emission-line region (BELR) "clouds" of the nuclear continuum radiation 
(e.g., ^ikora, Begelman, fc Rees 1994| ; piandford fc Levinson 1995|) or of the beamed jet 
emission itself ([Ghisellini fc Madau 19961) , as well as from the near-infrared emission of 
warm dust outside the BELR ( |Blazejowski et al. 20001) . As we discussed in §3.3| , if the 
emission from a radiative shock is dominated by IG scattering of external photons, then the 
synchrotron radiation component of the shock will have a relatively stronger contribution 
from regions of high values of B/B^q (see Fig. |Tn]).Q Most of the blazar emission models 
considered to date have, however, taken 7^ to lie below 7c (e.g., [Maraschi, Ghisellinij 
fc Gelotti 1992| ; Sikora et al. 1994; but see [Sikora et al. 20(JI| ) and thus correspond to 
nonradiative shocks. It is, nonetheless, likely that the cooling times could be quite short in 
these sources, as suggested, for example, by very rapid flares in blazars like 3G 279 (e.g., 
Wehrle et al. 1998|) , Mrk 421 (e.g., |Maraschi et al. 1999|) , and AO 0235+164 (e.g., [Kraus 
et al. 1999|) . Furthermore, there are indications that AGN jets may contain a significant, 
or even dominant, electron/positron component (e.g., [Reynolds et al. 1996^ [War die et al. 



1998[ ; ^ikora fc Madejski 2000[ ), which would naturally lead to high values of in the 



shocks. Given that the outfiows are likely magnetically driven (e.g., Blandford fc Znajek] 
1977[ ; [Blandford fc Payne 1981^ ), comparatively high values of may also be expected. 



It is thus conceivable that at least some of the synchrotron emission from these sources 
originates in radiative shocks in which the magnetic field is amplified by cooling-induced 
compression. 



^^The emission characteristics of GRB sources might also be affected by the interaction of the outflow 
with an external radiation field. For example, Lazatti et al. (2000) considered the possibility that the 
prompt 7-ray emission arises in the course of the Compton-drag deceleration of a relativistic outflow as it 
propagates through the debris of a progenitor massive star; in this scenario, magnetic fields play no role 
in the emission process. On larger scales, Madau, Blandford, & Rees (2000) proposed that backscattering 
of the GRB radiation by the comparatively dense massive-star environment could lead to a delayed MeV 
emission, again through the Compton-drag effect on the bulk flow. In the latter case, ERG emission by a 
"hot" (i.e., relativistic) electron/positron component of the flow could in principle also contribute to the 
observed radiation (cf. the analogous proposal for blazars by Ghisellini & Madau 1996). So far, however, 
there has been no discussion in the literature of the possible appearance of an ERG emission component in 
GRB sources. 
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5. Conclusions 

Gamma-ray bursts and their afterglows, as well as energetic flares from blazars and 
miniquasars, are commonly interpreted in terms of shocks in a relativistic outflow from a 
compact object that emit nonthermal (synchrotron and inverse Compton) radiation. In 
this paper we have investigated the structure of radiative shocks in such sources. A shock 
is "radiative" if the emission mechanism can tap into the bulk of the postshock internal 
energy and if the radiative cooling time behind the shock transition is shorter than the 
relevant dynamical time (the characteristic time for adiabatic energy losses). Such a shock 
is characterized by a radiative zone of flnite width where most of the energy dissipated in 
the shock transition front is radiated away. The structure of radiative shocks has been 
previously studied in the nonrelativistic regime in the context of the interstellar medium, 
where the dominant emission processes are thermal. One particularly interesting finding 
of these studies has been the strong effect that such a shock could have on the normal 
component of a magnetic field that is frozen into the gas, and, in turn, the potentially 
important role that the field could play in limiting the compression in the shock. When the 
shock is nonradiative, the compression (712/711) and resultant magnetic field amplification 
{B2/B1 = 77.2/^1) have a value of at most a few in the nonrelativistic regime (generalizing 
to a few times the shock Lorentz factor 71 in the relativistic case, with n and B measured 
in the fluid rest frame). Efficient cooling can lead to a strong enhancement of the density 
and the fleld amplitude in the radiative zone, but when the magnetic pressure attains 
equipartition with the thermal pressure, further compression is inhibited and subsequent 
cooling proceeds at nearly constant density. We have shown that, when synchrotron 
radiation dominates the cooling, the interplay between the cooling and the compression 
becomes more pronounced on account of the feedback effect between the field amplification 
and the emission process: a stronger magnetic field increases the emissivity, which in turn 
induces a larger compression that further amplifies the field. 

The cooling-induced compressional field amplification may have potentially significant 
consequences in sources that harbor radiative shocks: if the preshock magnetic field is low, 
this mechanism could in some cases increase it to a level where the synchrotron radiation 
becomes detectable, and if the preshock field is already moderately strong (so that it 
is within an order of magnitude of equipartition in the shocked but uncooled gas), this 
process could result in a considerable fraction 0.1) of the radiation being emitted from 
regions where B ?a B^q. When the shock is radiative but the field immediately behind 
the shock transition is below equipartition, the initial cooling would be dominated by the 
inverse- Compton process: either synchrotron self-Compton (SSC) or external-radiation 
Compton (ERC), depending on which seed- photon energy density is higher in the shock 
frame. In this case, a larger fraction of the synchrotron emission comes from regions with a 
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stronger magnetic field if ERC dominates SSC. 

For standard scenarios of GRB outflows and their environments, this fleld amphflcation 
mechanism may be relevant to internal shocks that produce the prompt 7-ray emission 
and possibly also to the reverse shock responsible for the optical "flash" and radio "flare" 
in sources like GRB 990123. It will likely not be of much relevance to the forward shock 
that gives rise to the GRB afterglow, although we pointed out that if the shock expands 
into a magnetized wind from a Wolf-Rayet star then no further field amplification would 
be required to account for the lowest B / Bcq values that have so far been inferred from 
afterglow observations. We also suggested that radiative shocks in which this mechanism 
operates could occur in blazars. 
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Fig. 1. — The speed (/52, solid line), proper speed {u2, dashed line), and magnetic-to-thermal 
pressure ratio (^2, dash-dotted line) immediately behind the shock transition, as a function 
of the proper Alfven speed uai- The asymptotic values for uai much larger or much smaller 
than 1 are given in the figure. In order for a shock wave to exist, uai must be smaller than 
the upstream speed ui. 
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Fig. 2. — The asymptotic proper speed Umm (in units of U2) as a function of rj for several 
values of uai- 
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Fig. 3. — The asymptotic proper speed Wmin (in units of ^2) as a function of uai for several 
values of 77. 
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Fig. 4. — The proper speed -u, normalized by the postshock value ^2, as a function of the 
normalized distance ^ = x/xc for 1] = 0.9, 6 = 0, and several values of uai {solid lines). The 
dashed lines represent the corresponding asymptotic values (Mmin). The boldface dash-dotted 
line depicts the limit uai 0, and the boldface dotted line is the corresponding asymptotic 
value. An asterisk is placed on the sohd lines and on the dash-dotted line at the point where 
u has been reduced to 90% of its asymptotic value: {u — Uiam)/iu2 — Umin) =0.1. 
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Fig. 5. — The proper speed u, normalized by the postshock value M2, as a function of the 
normalized distance ^ = x/xc for uai = 0.05, 6 = 0, and several values of r] {solid lines). 
The dashed lines and the asterisk symbols have the same meaning as in Fig. ^. 
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Fig. 6. — The pressures of the electron/positron (pe) and proton {pp) components, the total 
thermal pressure (p = Pe + Pp), the magnetic field pressure {pb — B'^/8n) and the total 
pressure (ptot — P + Pb), all normalized by the postshock thermal pressure p2, as a function 
of the normalized proper speed u/u2 for 77 = 0.9 and uai — 0.2. 
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Fig. 7. — The minimum value of rj for which the magnetic field in the radiative zone reaches 
equipartition with the total thermal pressure, as a function of uai- For uai > Ucq the 
magnetic field amplitude just behind the shock transition is already above equipartition. 




Fig. 8. — The minimum value of rj for which the magnetic field in the radiative zone reaches 
equipartition with the total thermal pressure, as a function of 82, the ratio of the postshock 
magnetic pressure to its local equipartition value. 



-43- 




Fig. 9. — The fraction of the postshock internal energy still left in the electron/positron 
component as a function of the postshock magnetic-to-thermal pressure ratio, ^2: (i) at 
the point where the magnetic field attains equipartition with the thermal pressure, S — 1 
{solid curves), and (ii) at the point where the magnetic pressure is 1% of the thermal 
pressure, S = 0.01 {dash- dotted curves). The curves in each of the two cases correspond 
to 77 = 0.5, 0.75, 0.9, 0.95, and 1, respectively, from right to left. 
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Fig. 10. — The solid lines show the fraction of the synchrotron emission that is emitted in 
a region where the magnetic field is larger than a given value, Psyn{> -B)/-Psyn (total), as a 
function of the magnetic field strength (normalized by its equipartition value). The dashed 
lines show the fraction of the total emission that is contributed by the synchrotron process, 
Psyn(> -B)/-Ptotai(> B). The results are shown for three sets of parameters: (A) 6 = 0, 77 = 1, 
UAi = 0.03; (B) 6 = 0, ?7 = 1, uai = 1; (C) b = 100, 77 = 1, uai = 0.03. 



